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Abstract. Interpreting the J/ip suppression reported in nucleus-nucleus collisions at SPS and RHIC re- 
quires the quantitative understanding of cold nuclear matter effects, such as the inelastic rescattering of 
J/ip states in nuclei or the nuclear modification of parton densities. With respect to our former Glauber 
analysis, we include in the present work the new PHENIX d-Au measurements, and analyze as well all 
existing data using the EPS08 nuclear parton densities recently released. The largest suppression reported 
in the new PHENIX analysis leads in turn to an increase of a J/ifiN from 3.5 ± 0.3 mb to 5.4 ± 2.5 mb using 

proton PDF. The stronger x-dependence of the G A /G P ratio in EPS08 as compared to e.g. EKS98 shifts 
the cross section towards larger values at fixed target energies (x2 ~ 0.1) while decreasing somehow the 
value extracted at RHIC (x2 ~ 10 -2 ). 

PACS. PACS-key discribing text of that key - PACS-key discribing text of that key 



1 Introduction 



The suppression of heavy-quark bound states in heavy-ion 
collisions due to Debye screening is known to be a sensitive 
probe for quark-gluon plasma (QGP) formation [lj. How- 
ever, reactions involving heavy nuclei introduce "cold" nu- 
clear effects which are not due to QGP formation but that 
affect J/ip production nonetheless. Among them, the nu- 
clear modification of the parton densities (nPDF) may 
play a role in the nuclear dependence of J/ip production. 
Another effect is the inelastic rescattering of the J/tp state 
in a nuclear matter, the so-called nuclear absorption. It is 
therefore crucial to have a quantitative understanding of 
these cold nuclear effects in order to get a quantitative 
understanding of the J/ip suppression reported in heavy 
systems at SPS an d RHIC [4] and, therefore, a reli- 
able to interpretation of the observed suppression. 

Nuclear absorption is expected to be the dominant 
source of J/ip suppression not only in peripheral heavy 
nucleus-nucleus (AA) collisions but also in hadron(photon)- 
nucleus reactions, which are dominantly sensitive to cold 
nuclear effects. Experimentally, a large variety of J/ip pro- 
duction off nuclear targets have been measured at various 
colliding energies (SPS |5|lrJ|T7j. FNAL [5|ig|HD|UT] . HERA- 
B [12], RHIC [13], SLAC [H] and NMC [15]). A global 
analysis of all available measurements of J/ip production 
in nuclear target allows for a study of nuclear absorption 
effects and to quantify the strength of that mechanism, 
monitored essentially by one physical parameter, the J/ip- 
nucleon inelastic cross section, a , , , M . In the following, the 
extraction of ct j/VjN assuming nuclear parton modifications 



is presented within the framework used in a previous anal- 
ysis [16] . 

In this paper, the most recent measurements performed 
by the PHENIX experiment in d-Au collisions [17] are 
analyzed. Moreover, additional results using the EPS08 
nPDF parametrization [18j are given. 



2 Extracting a^ N 

This section gives a brief description of the method fol- 
lowed in this analysis: the model used to describe the data 
selected, the nuclear parton distribution implementation, 
the data sets and finally the fitting method. A more de- 
tailed description of the method can be found in [16] . 



2.1 Model 

The various J/ip production channels in the different re- 
actions in the data sample are the following: 



,7*)+ J/ip + X 



(1) 



The J/ip production cross section crP™^ in hadronic colli- 
sions is determined within the Colour Evaporation Model [19] 
(CEM) at leading order (LO). The PDF in the hadron pro- 
jectiles are taken from the LO parametrization SMRS for 
the pion [20] and CTEQ6L for the (anti)proton [21]. Since 
only cross section ratios are analyzed in the following, the 
results from this analysis are almost independent on the 
specific choice of the proton PDF parametrization. 
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The survival probability S a bs (A CT j/^ N ) 01 J/ip states 
propagating in a nucleus A - i.e. the probability for no 
inelastic interaction - is given in the Glauber model by [H] 



^abs (A, (J j/w, N ) — 



1 



(A 1) a 3 



db 1 



e -(l-lM) T A (b) ^ N)(2) 



with the thickness function T A (b) 



r A (b) 



d z p(b, z). 



(3) 



It depends on the atomic mass number A of the nu- 
cleus and the J/ip—N inelastic cross-section, <J J/4>N - The 
observed J/ip production as a function of the longitudinal 
momentum fraction x, then is: 



dcr 



dx 



S'abs (A , CTj/v.n) x 



d(T prod 



dx 



(4) 



In this current analysis, the cross section ratios R th of 
heavy (A) to light (B) nuclei are considered: 



Sda(/»,7*A -> J/ipX)/dx 
Ada^fB J/ipX)/dx 



(5) 



Note that since only ratios of cross sections at the same 
energy are used, most uncertainties regarding the J/ ^> pro- 
duction cross sections cancel. 

It is worthwhile to note that formation-time effects are 
neglected, in the sense that the question of which state 
actually propagates through the nuclear matter is not ad- 
dressed. Also, the feed down from higher mass resonances 
is not taken into account. Consequently, c //VjN has to be 
seen as an effective absorption parameter resulting from 
an average of the cc and J/ip, Xc and ?p' interaction with 
nucleons, rather than the genuine J/tp-N inelastic cross 
section. 



2.1.1 Nuclear parton distributions 

Partons in bound nucleons show noticeably different mo- 
mentum distributions as compared to those in free pro- 
tons. This modification is quantified by R(x, Q 2 , A) as a 
function of Bjorken variable, the square of the momen- 
tum transfer Q 2 and the nucleus size A in the following 
formula : 



R l (x,Q 2 ,A) = f?(x,Q 2 ) I Aff(x,Q 2 ) 



(6) 



where fi and /• describe respectively the distribution of 
parton i in a proton and in a nucleus. 

Since J / ip is predominantly produced via gluon fusion 
inp-A collisions)]] at high energy (y/s^A >20 GeV) its pro- 
duction is affected by the modification of the gluon distri- 
bution in nuclei. Several DGLAP analyses 23,24,25 aim 



1 In 7r ± -A and p-A collisions, the scattering of a valence 
antiquark from the projectile with a valence quark from the 
target is favoured. 




Fig. 1. The ratio of the gluon distribution in a gold nucleus 
over that in a proton, R^ u (x, Mji^), plotted as a function of 
the Bjorken variable using the nDS, nDSg [23], EKS98 [55], 
HKN [24] and EPS08 [18] parametrizations. The bands in- 
dicate the typical rr-range probed by J /if) production in the 
NMC, SPS, FNAL, HERA-B, and RHIC experiments (top to 
bottom). 



at the extraction of the ratios Ri(x, Q 2 , A) from DIS and 
Drell-Yan data. However, given the indirect constraints 
in the gluon sector (through scaling violations), the ratio 
R g is poorly determined. Figure [1] shows the gluon dis- 
tribution in a Au nucleus with various parametrizations 
available as a function of x. 

The shaded band area shows the kinematic range of 
the J/ip production (at LO) for various experiments, NMC 
(green), SPS (blue), FNAL (red), HERA-B (orange) and 
RHIC (purple), from top to bottom. One can observe that 
J /if) production is affected by mainly two effects, the anti- 
shadowing (i? Au > 1 at 2-5 x 1CT 2 < x < 0.3) for SPS, 
FNAL and HERA-B and the shadowing effect (i? Au < 1 
at x < 10~ 2 ) at RHIC. A strong anti-shadowing effect 
increases J/tp production in nuclei with respect to the 
(binary scaled) production in p-p collisions, leading to a 
cross-section ratio larger than 1. This enhanced produc- 
tion will in turn be compensated by an increase of the 
fitted nuclear absorption cross section (SPS, FNAL and 
HERA-B). Conversely, a strong shadowing effect tends to 
deplete the nuclear absorption cross section (RHIC en- 
ergy). 

In this work, the EPS08 parametrization (magenta, 
dotted-dashed-dashed line) is added in the analysis. This 
nPDF set exhibits a strong anti-shadowing effect in com- 
parison with the previous EKS98 distributions, and the 
anti-shadowing region is also slightly shifted to higher x 
values (2 x 10~ 2 < x < 0.3). In addition, the shadowing 
effect is much stronger than in EKS98, due to the inclu- 
sion in the analysis of these authors of the recent RHIC 
data. 



2.2 Data sets 

Since factorization between J/ip production and the sub- 
sequent inelastic interaction is assumed in the present 
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analysis, both hadroproduction (using pion, proton, an- 
tiproton and deuterium beams) and leptoproduction data 
are analyzed. 

The detailed data selection list can be found in [16] . 
Concerning hadroproduction measurements, the projec- 
tiles used were mainly protons (NA3, NA38, NA50, E866, 
HERA-B, PHENIX), but also pions (E537, NA3, E672), 
antiprotons (E537), and deuterium nuclei (PHENIX). The 
range of colliding energy is ^/snn = 15-200 GeV. 

As mentioned previously, PHENIX data have been re- 
analyzed [17j , the measurements in <i-Au are now normal- 
ized with respect to higher-statistics p-p collisions mea- 
surements performed at positive/negative rapidity instead 
of an average of measurements at positive/negative rapid- 
ity. The i? exp ratios are now smaller than in the previous 
analysis [4] for each rapidity region. Concerning the uncer- 
tainties, the precise p-p measurements also lead to slightly 
smaller statistical uncertainties. However, since the data 
are not taken during the same year and with the same 
configuration, the systematic uncertainties which used to 
cancel in the ratio are now larger. 

As for leptoproduction experiments, the NMC data [15] 
are selected. The virtual-photon energy v ranges from 40 
to 240 GeV in the laboratory frame, corresponding to 
7*-N centre-of-mass energies ^fs = 8-21 GeV and the 
Bjorken-x range probed in the gluon distributions of the 
nuclear targets is x = 0.05-0.15 in the NMC kinematics. 

After the data selection, the i? Gxp ratios of heavy (A) 
to light (B) nuclei are determined and compared to the 
i? th ratios. In order to avoid too large systematic errors in 
the experimental ratio, both reactions on targets A i and 
B are required to be taken from the same experiment and 
at the same centre-of-mass energy. For each experiment, 
the uncertainties on the i? exp ratios are then separated as 
follows: 



through: 



R™ p ±a i ±0R e i 



(7) 



where Cj represents the uncorrelated errors (statistical and 
the uncorrelated systematic errors, added in quadrature), 
and j3 corresponds to the normalization correlated error, 
often coming from the fact that cross sections in different 
nuclei are all normalized to the same light target (hence 
with an error common to all Ri). 



2.3 Fitting method 

The J/tp-N inelastic cross section is extracted, for each 
experimental sample £ with n e data points, from the min- 
imization of the x 2 , function [26] : 



1=1 



R 



cxp 



~ R i (°"j/*n) 



V 2 /M, (8) 



computed from Eq. l[5]), depends explicitly on the free, but 
positive, parameter, cTj,, N . The correlated normalization 
error /?, on the data point i enters the V and M in Eq. © 



i=l 



M = 1 +/3 



2 \ - 



The la error Sa j/j/jN on the fitted parameter a J/<jiN is de- 
fined so as to increase % 2 by one unit from its minimum: 

A X 2 ee ± 5a JN J - X 2 _ = 1. (9) 



3 Determination of (x . from each 
experiment 

The analysis using new PHENIX results are compared to 
the previous analysis in Table [1] Since the i? exp ratios in 
the PHENIX new analysis [T7] are smaller than in [1] , the 
nuclear absorption cross section obtained in this work is 
now higher than previously, by roughly 2 mb. The cross 
sections now vary from 2.5 ±2.2 mb (with nDSg) to 5.4 ± 
2.5 mb (proton PDF) with various nPDF parametriza- 
tions. 

The extracted nuclear absorption cross sections us- 
ing the EPS08 nPDF paramctrization arc shown in Ta- 
ble[2]for all individual experiments. These results are com- 
pared with the results obtained previously in using 
the EKS98 parametrization. Because of the more pro- 
nounced shadowing in the EPS08 parametrization than 
in EKS98, the extracted nuclear absorption cross section 
at RHIC energy decreases by ~ 30%. The anti-shadowing 
effect is also more pronounced in the EPS08 parametriza- 
tion, leading to an increase by roughly 10-20% of the nu- 
clear absorption cross section at the energies of the SPS, 
FNAL and HERA-B experiments. Finally, note the sig- 
nificant increase (+50%) from EKS98 to EPS08 using the 
E866/Nusea data samples. However the large x 2 /ndf ~ 20 
obtained is large, as stressed in [2"7] . 

The nuclear absorption cross section for PHENIX de- 
pends on the strength of shadowing in each nPDF set: 
a strong shadowing parametrization leads to a decrease 
of the nuclear absorption cross section. On the contrary, 
when using a proton PDF or a nPDF with no (or weak) 
shadowing effect, the nuclear absorption cross section is 
then higher to compensate for the (weak) suppression due 
to gluon shadowing. 

Since the energy in the J/^-nucleon (or cc-nucleon) 



system, given by ^fsjj 



ij/^/y/x^, is directly related 



J/jpN — J/V 

to the momentum fraction X2, one could expect the ex- 
tracted o~j,, N cross section to be a scaling function of xi . 
As discussed in [16], there is no real ^-dependence ob- 
served within this framework using a proton distribution, 



The result obtained using the E672 measurements does not 
exhibit this increase, the x region probed (x ~0.03) lies in the 
anti-shadowing region in EKS parametrization while within the 
EPS parametrization, this x region is at the boundary of the 
shadowing region with a weak effect of the nuclear distribution. 
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Table 1. The J/tp-'N cross section extracted from the new re-analyzed PHENIX data versus previous analysis using the proton 
and various nuclear parton density parametrizations. The x 2 / n df and the % 2 probability are also shown. 





Previous results using 13 
ct j/^n ( mb ) X 2 /ndf Probability 


New analysis using [17 
a j/i,N ( mb ) X 2 /ndf Probability 


Absolute change 


proton 


3.5 ± 3.0 


1.7 


0.79 


5.4 ± 2.5 


0.84 


0.93 


+1.! 


) mb 


nDS 


3.1 ± 2.6 


1.4 


0.84 


5.1 ± 2.5 


0.69 


0.95 


+2.0 mb 


nDSg 


0.6 ± 1.9 


0.8 


0.93 


2.5 ± 2.2 


0.27 


0.99 


+ 1.S 


) mb 


HKN 


1.5 ± 2.3 


1.3 


0.86 


3.2 ± 2.3 


0.56 


0.97 


+1.7 mb 


EKS98 


1.3 ± 2.0 


0.6 


0.93 


3.1 ± 2.2 


0.12 


1.00 


+U 


1 mb 


EPS08 


1.3 ± 2.5 


1.5 


0.83 


2.2 ± 2.2 


0.37 


0.98 


+0.! 


) mb 



Table 2. The J/V'-N inelastic cross section, x 2 /ndf extracted from each data sample using EKS98 and EPS08 parametrizations 
for the nuclear PDFs. 



Exp. 


<r??f M (mb) 


XEKs/ ndf 


a*™ (mb) 


X E Ps/ n df 


Relative change 


E537 


8.2 ± 1.1 


1.9 


9.0 ± 1.2 


1.86 


+10% 


NA3 


4.6 ± 0.2 


1.2 


5.2 ± 0.2 


1.32 


+13% 


NA38 


7.9 ± 0.8 


3.2 


9.0 ± 0.8 


3.07 


+14% 


NA50 


6.8 ± 0.5 


0.3 


7.8 ± 0.5 


0.31 


+15% 


E672 


11.6 ± 6.3 


0.6 


10.0 ± 5.8 


0.61 


-14% 


E866 


5.3 ± 1.7 


6.5 


8.0 ± 3.7 


20.4 


+51% 


HERA-B 


4.2 ± 1.5 


0.9 


5.1 ± 1.5 


0.8 


+21% 


PHENIX 


3.1 ± 2.2 


0.12 


2.2 ± 2.2 


0.37 


-29% 


NMC 


< 1.6 


0.5 


< 2.00 


0.35 


+25% 















or using nDS, nDSg, EKS98 and HKN nuclear parton dis- 
tribution. For completeness, Figure [5] shows the nuclear 
absorption cross section ct //i/jN as a function of X2 using 
the EPS08 nPDF. In the region of s 2 ~ 0.1, one can 
observe that the spread of extracted a J/l/jN reported us- 
ing the other nPDF sets persists. Interestingly, it also ap- 
pears that using EPS08 leads to some decrease of a Jf 
from fixed-target to RHIC energies, indicating possible 
formation-time effects at small x-i. Also, in the previous 
analysis [16] , a similar trend has been observed when us- 
ing EKS08 paramctrization. However, the error bars are 
too large to make any firm conclusion for EKS08/EPS08 
and other nPDF sets. Note that higher-twist production 
mechanisms may very well have a different kinematic de- 
pendence; this is for intance for the case for the intrinsic 
charm model which naturally exhibits a Feynman-x scal- 
ing (see e.g. (29)- However, we expect its contribution to 
be small at low \xf\ which we consider here. 



4 Global fit and discussions 

In the following, a global fit is performed assuming that 
the (Jj.^jj dependence on energy is weak. The detailed 
method for the global fit is described in [16]. The la error 
is rescaled, 



da 



S x 5a 



7/-0N ' 



where the factor S is defined by: 

^— if x7ridf>l, 
n — 1 



S = 



(10) 

(11) 



Q 

E 



12 



10 



• E537 

T NA3 

★ NA38 
O NA50 
A E672 
■ E866 

▲ HERA-B 

* PHENIX 
A NMC 



EPS08 




10 



10 



10 



Fig. 2. The J/t/j-N cross section extracted from each data set, 



using EPS08 {a 



nuclear parton densities as a function of 



X'Z- 



with n data points and S = 1 otherwise. The J/ip-N cross 
section is systematically determined from the individual 
data samples. 



The extracted a 



is then determined from the mini- 
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mization of the weighted x 2 function: 

N 

x 2 K* N ) = £ S ^ X?K/*n)- (12) 

1=1 

with the individual \e f° r each experimental data sam- 
ple. This global fit analysis will thus favour data sets 
with a small individual x 2 /ndf. The results obtained from 
this global fit using a proton PDF and various nPDF 
parametrizations (nDS, nDSg, HKN, EKS98 and EPS08) 
are summarized in Table [3] These results include the re- 
cent PHENIX results [T7] already mentioned. The x 2 /ndf 
from these fits varies from 1.4 to 1.7. 

The spread of Cjf^f quoted in Table [3] directly reflects 
the present lack of knowledge of the (gluon) nuclear den- 
sities. Taking the nDS parametrization as the default set, 
the cross section extracted in this analysis is: 

<r J/ ^ N = 3.5 ± 0.2 (stat.) ± 2.6 (syst.) mb, (13) 

where the systematic error quoted here only comes from 
the uncertainties of the nPDFs. Clearly, a better deter- 
mination of u nPDF could only be achieved when these are 
more tightly constrained by experimental data. 

Figure shows the fitted Cj,^ n from this work com- 
pared to another global analysis (Gerschel and Huffier 
(GH) in [29], Kharzeev et al. (KLNS) in Ref. [30]) as well 
as with the extracted <t j/i(jN by NA50 [7] and PHENIX [17] 
from their data. Both results by Gerschel and Hiifner (GH) 
in [29] and Kharzeev et al. (KLNS) in Ref. [30] are sig- 
nificantly higher than the o^. cross section presented in 
this work. These differences are believed to be mainly due 
to the different data sets used in the global analysis. A 
more detailed discussion can be found in [16j . 



the collaboration also published the value of cr^. us- 
ing the nDSg and EKS98 parametrizations. The results 
are compatible with results presented in this work within 
the error bars. When using the EKS98 parametrization, 
PHENIX results are 2.%t\ A vs. 3.1 ± 2.2 mb in this work 
and when using nDSg, PHENIX results are 2. 2^'!? mb vs. 
2.5±2.2mb. 

5 Summary 

In this work, a re-analysis of the nuclear absorption cross 
section using the new PHENIX results within the frame- 
work described in [16j is presented. The largest suppres- 
sion reported in the new PHENIX analysis leads to an 
increase of a ... from 3.5 ± 0.3 mb to 5.4 ± 2.5 mb using 
the PDF of the proton. The cr^. obtained in this work 
is also compatible within the uncertainties with the value 
determined from the PHENIX analysis of their measure- 
ments. It is worthwhile to note that RHIC has provided 
high-statistics d-Au. collisions during the 2008 year data 
taking, and the analysis of this new set of data should 
allow for a more precise measurements of the R cxp ratios. 

In addition, an analysis of the c 7/i/jN nuclear absorp- 
tion cross section is performed using the EPS08 nPDF set; 
it is presented for each individual experiment. The strong 
shadowing and anti-shadowing effects described by this 
parametrization induce in turn a possible xi dependence 
of the cr J/</)N cross section, leading to a smaller nuclear 
absorption at low xi (RHIC energy) and increasing it at 
large x 2 (SPS, FNAL, HERA-B). However, the discrep- 
ancy of the extracted a - N observed at large x (x ~0.1) 
persists. Finally, a global fit including the new PHENIX 
results and the EPS08 parametrization is presented in this 
work. 
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Fig. 3. The J/i/)-N cross section extracted in the global fit in 
this work compared to previous analyses. 

Compared to NA50 [7] analysis, the cr N extracted 
on their measurement is compatible with the individual 
a J/iiN extracted from this work using a proton parton dis- 
tribution, namely er NA50 = 4.2±0.5 mb vs. 4.7±0.5 nib in 
this work. In the new analysis of the PHENIX data [17], 
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